Time-dependent stabilization in AdS/CFT 



Roberto Auzzi, Shmuel Elitzur, Sven Bjarke Gudnason and Eliezer Rabinovici 

Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 
E-mail: auzzi (at) phys .huj i . ac . il, elitzur (at) vms .huj i . ac . il, 
gudnason(at)phys .huj i . ac . il, eliezer (at) vms .huj i . ac . il 

Abstract: 

We consider theories with time-dependent Hamiltonians which alternate between being 
bounded and unbounded from below. For appropriate frequencies dynamical stabilization 
can occur rendering the effective potential of the system stable. We first study a free field 
theory on a torus with a time-dependent mass term, finding that the stability regions are 
described in terms of the phase diagram of the Mathieu equation. Using number theory we 
have found a compactification scheme such as to avoid resonances for all momentum modes 
in the theory. We further consider the gravity dual of a conformal field theory on a sphere 
in three spacetime dimensions, deformed by a doubletrace operator. The gravity dual of 
the theory with a constant unbounded potential develops big crunch singularities; we study 
when such singularities can be cured by dynamical stabilization. We numerically solve the 
Einstein-scalar equations of motion in the case of a time-dependent doubletrace deformation 
and find that for sufficiently high frequencies the theory is dynamically stabilized and big 
crunches get screened by black hole horizons. 
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1 Introduction 



One of the main challenges facing string theory is to uncover what is to actually become 
of various types of classical gravitational singularities. The AdS/CFT framework is a 
particularly useful setup for such a study. On the one hand it provides complete non- 
perturbative information and on the other hand, decays from a false vacuum in AdS [1] 
generically result classically in a big crunch. This was studied in [2-13]; in particular it was 
argued in [10, 13] that there are cases were the singular nature of a crunch is a question of 
observables and conformal frame used to describe them. The description of a big crunch as 
a boundary theory can be singular in one frame and non-singular in another. In the bulk the 
singularity is not resolved and yet infinite entropy is encoded on the boundary. In one frame 
the Hamiltonian is always bounded and in the other it is bounded as long as it exists. There 
are other more radical crunches which involve infinite energies and unbounded Hamiltonians 
on the boundary [2-5, 10, 13]. In those cases there is no clear expectation that the crunch 
would or should be cured. In this work we analyze boundary and bulk behavior in a class 
of cases in which the boundary Hamiltonian's properties are somewhat midway between 
the bounded and unbounded cases. We study here a class of time-dependent Hamiltonians 
which are unbounded "half of an infinite time and bounded during the other "half of 
time. This in an oscillatory manner. A butterfly flapping its wings. Such systems are 
interesting to discuss in their own right. In some range of parameters it turns out that the 
energy of such systems is effectively bounded from below. We will survey how the stability 
or instability features of these boundary potentials manifest themselves in the bulk. When 
crunches would occur, would black holes form to shield the singularities? In particular, one 
could pose the question if it is possible to turn the mentioned big crunch into a near-crunch 
by means of such a dynamical stabilization. As a prototype example we will consider in 
this paper the same unbounded potential, however, multiplied by an oscillating function 
such that half of the time the potential is bounded and half it is not. Concretely we will 
consider an unbounded potential of the dual boundary theory, oscillating with frequency 
oj as follows 



where (j> is the field in the problem at hand. An example of dynamical stabilization is 
Kapitza's pendulum [14] which is stabilized for frequencies above some critical frequency 
uj > lj c . This was done using the method of separation of time scales. 

In order to acquire some intuition we will first consider the example of a free field 
theory. For the free field theory it will prove convenient to break down the problem to 
the zero mode and higher momentum modes. This leads us to study analogous quantum 
mechanics problems before considering the field theory cases at hand. We find in the 
case of a free massive field theory, when compactified on a torus, that it will indeed be 
stabilized when the (angular) frequency is larger than an order one number times the 
inverse compactification radius. There is though an unwanted side effect coming about by 
means of our treatment, i.e. dynamically stabilizing the system with oscillatory behavior, 
will in general lead to a secondary effect of resonances. This is a different type of instability 
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not present in the initially static unbounded potential. An analysis in terms of number 
theory shows that in the case of the free field theory it is possible to avoid the resonances 
by choosing certain values of the compactification radius. 

The free massive field theory has the Mathieu equation as its equation of motion, which 
shows up many places in time-dependent physical problems. Tacoma bridge and the Paul 
trap are classic examples. Other more recent phenomena involving the Mathieu equation 
include reheating in the context of cosmology [15-18], relativistic ion traps [19] as well as 
trapping of particles using lasers [20] . 

In the interacting field theory the story is somewhat more complicated and we do not 
have concrete proofs at hand. Intuition and the method of separation of time scales provide 
evidence for the same to occur in this case as well. A difference between the free and the 
interacting theory is the expected thermalization in the latter. However for sufficiently high 
frequencies we have indications from the zero mode by means of the quantum mechanics 
problem as well as from the bulk that the boundary theory is dynamically stabilized. 

Concretely we solve the Einstein-scalar equations numerically in global AdS4 deformed 
by a time-dependent doubletrace operator finding that the big crunch may or may not 
appear depending on the choice of frequency. In all cases we start by having a smooth 
initial field configuration given by a Coleman-de Luccia (CdL) instanton [1] (restricted to 
the equator). 

Most of the numerical studies in AdS space concerning Einstein-scalar equations so far 
have focused on a massless scalar field (e.g. [21-25]). In order to be able to introduce a 
(non-irrelevant) multitrace deformation [26-28] in the case of AdS^+i, it is required that 
the mass squared of the scalar obeys 

d 2 2 d 2 
__< m 2 < __ + 1) (L2) 

such that both possible fall offs of the scalar field 

^7a7 + 7a7' A ± = -±^/ T + m2, (1.3) 

are normalizable [29]. The multitrace deformation is then realized as a boundary condition 
relating a and /3. Hence it proves crucial to consider an AdS field with a mass squared in the 
range (1.2) in order to allow the introduction of an external time-dependent potential. A 
doubletrace deformation in the boundary theory (which is proportional to H (t)0 2 /2, where 
O is an operator dual to the bulk field 4>) corresponds to the boundary condition (3 = H(t)a. 
We numerically solve the Einstein-scalar equations subject to the latter boundary condition 
with H(t) = /locos(wt). This enables us to monitor the circumstances under which black 
holes are formed. We are particularly interested in the oj dependence of the formation. 

Having set the scene we will start out by recalling the method of separation of time 
scales classically (Sec. 2) and quantum mechanically (Sec. 3) which we generalize to field 
theory in Sec. 4 for both a free massive field theory as well as for a quartically interacting 
theory. In Sec. 5 we study in detail the bulk manifestations of the stabilities and instabilities 
on the boundary. We will present evidence for big crunches which in some cases are indeed 
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shielded by black hole horizons. This provides evidence that thermalization indeed occurs 
and is responsible for the dynamical stabilization in the case of interest. Specifically we 
will consider a system in global AdS4 being dual to a 4D CFT on a sphere deformed by a 
doubletrace operator. 

2 Classical mechanics 

In this section we will review the merits of the method of separation of time scales. This 
method was used by Kapitza [14] to understand the dynamical stability of a vertically 
rapidly oscillating suspension holding a pendulum. The system is stabilized with the pen- 
dulum pointing straight upwards due to the kinetic energy of the very fast oscillations. 
Whenever the pendulum is slightly displaced from the point of equilibrium the effective 
potential pushes the pendulum back into the straight upwards position with damped slow 
oscillations. 

Here we review the derivation of the effective Hamiltonian (viz. the Hamiltonian de- 
scribing the drift part of the system) of a rapidly oscillating potential using the method of 
separation of time scales following [30-33] but specialized to the Hamiltonian 

2 

H = — - V (x) cosfwt) . (2.1) 
2m 

We decompose the particle x(t) into a slowly moving part X and a rapidly moving part £ 
as 

oo 

x(t)=X(t)+£(X,X,u)t), Z{X,XM) = ^T J u~ i Zi{X,XM), (2.2) 

i=i 

where £ is periodic with vanishing average in oot and the functions £j are chosen such that 
X is independent of oot. From the equation of motion 



m 



[X+'i) = V '(X + £)cos(u,i), (2.3) 



we can use the chain rule on derivatives of £ and in turn determine £j order by order in 
1/uj. Doing so up till order a; -4 , it is possible to determine the drifting part of the equation 
of motion, i.e. the terms that do not average to zero and hence are not absorbed into the 
£jS. Expanding out the left hand side of (2.3) does not give any non-vanishing terms upon 
time averaging (by definition), hence we need to expand the right hand side 

Vq(X) cos(wt) + £Vq(X) cos(wt) + -£V (3) (X) cos(wi) + • • • (2.4) 

The non-vanishing terms to 4th order upon time averaging are 

^ W'PO cosM + ^UV^(X) cosM) , (2.5) 

where the bar denotes average with respect to time: x = d(u)t)x / (2iv) . Calculating the 
coefficients of £ 

^cos(^) 2V£Xsm(ut) V{V» C0s(2u>t) 3F (3) X 2 cosjut) 5 

S — 5 1 S 1 3 T~A 1 4 T U V W ) ) ( Z -V 
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we obtain the drifting part of the equation of motion 

mjt = "is? W + 2^'^ 2 + °<^> ■ (2 - 7) 

This equation of motion can be obtained from the following effective Hamiltonian 



with P being the momentum conjugate of X. Physically we can explain this effective 
Hamiltonian as follows. Since our choice of potential averages out to zero in time, to 
leading order the theory is free. The correction to the potential equals |m£ 2 up to 4th 
order in 1/oj, i.e. it is the kinetic energy of the rapid oscillations [30] 1 . This leading order 
potential is always confining. Notice that the dynamics of the drift degree of freedom X 
happens to be conservative in spite of the fact that the original system has an explicit 
time-dependent Hamiltonian (2.1). For an alternative derivation of the effective potential 
and an explanation of the effective conservation of energy, see app. C. 

2.1 Parameter counting 

Considering a monomial of the coordinate in a one dimensional classical mechanics system 

2 

H=%-- X n+1 x n cos(wt) , (2.9) 
2m 

we can infer that A has dimension of mass. Performing a rescaling we can write 

dx 1 A\ n+1 L n ~ 2 



H = A 



v 2 

V --q~x n cos(2r) 



V 



— , A = -mu) 2 L 2 



dr ' 4 muj 2 

(2.10) 

where A (with unit of mass) sets the overall scale of the problem and q is a dimensionless 
parameter determining the rapidness of the oscillations (small q corresponds to fast oscil- 
lations). Finally x = x/L and r = ut/2 are dimensionless spatial and time coordinates, 
respectively. This shows that the system at hand has two parameters to dial (for fixed 
boundary conditions). It is not possible to use A as an expansion parameter while it is 
indeed possible to expand in small q — it is exactly the 1 /to expansion. 

2.2 Quadratic potential 

Let us consider the n = 2 case of the Hamiltonian (2.10), i.e. the case of the quadratic 
potential. Applying the general expansion, the effective Hamiltonian using eq. (2.8) for 
Vb = X 3 x 2 reads 

H cS = f- + ^- 2 X 2 + ^jP 2 + 0(^) . (2.11) 

The last term can be interpreted as modifying the effective mass with respect to the bare 
mass due to the oscillations. The potential under consideration is very special though. In 



1 At the 4th order in 1/oj, \m£, 2 gives the effective potential only up to a total derivative. 
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fact it is not needed to make any expansion as the classical problem is completely solvable. 
The equation of motion coming from eq. (2.10) with n = 2 (dropping the tilde) reads 

X - 2qcos(2r)x = 0, (2.12) 

which is the Mathieu equation (4.2) with a = 0. Let us dwell on this classical problem for 
a while. The general solution to this linear ordinary differential equation is given in terms 
of the so-called Floquet solution, which can be expressed as 

F( T ) = e iuT P{T) , (2.13) 

where P(r) is a periodic function with period tt and v is the characteristic exponent which 
is in general complex. Hence, the existence of an imaginary part of the characteristic 
exponent implies that the solution either blows up in the infinite future or the infinite past. 
The general solution is given by a linear combination of F{r) and F{— r) and therefore, 
unless one of the coefficients is tuned to zero, the solution will blow up in the future if the 
imaginary part of v is non-zero. 

Now we will first pose the question of stability. Let us define stability by boundedness 
of the generalized coordinate x, i.e. one could interpret it as the classical "particle" not 
running off to infinity exponentially fast. The answer in the classical case is given by the 
phase diagram of the Mathieu equation [34], see fig. 1. The first few stable bands are 
shown in table 1. As compatible with physical intuition, the stable region is roughly only 



band # s 


Qs, start 


Qs, end 


1 





0.908046 


2 


7.51361 


7.57969 


3 


21.2986 


21.3032 


4 


41.9272 


41.9274 


5 


69.4284 


69.4284 



Table 1. Stable bands in the variable q of the Mathieu equation (4.2) for a = 0. 



given by small q < 0.9, except for very narrow bands in q. According to eq. (2.10) large 
frequencies correspond to small q. Thus when the potential is oscillating very rapidly, the 
forces practically cancel and we are left with a "free particle". This agrees with intuition 
as well as the 1/co expansion. For an estimate of the amplitude of the particle position, see 
app. B. 

2.3 Quartic potential 

Turning to the n = 4 case of the Hamiltonian (2.10) and applying the expansion we obtain 
for Vq = A 5 x 4 the effective Hamiltonian 

H cS = f - + + ^V^P 2 + 0{^) . (2.14) 

To our knowledge this case is not exactly solvable and hence we have nothing to add to 
the results of the 1/uj expansion, i.e. the effective potential is confining. 
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Figure 1. Phase diagram of the Mathieu equation (4.2). The stable regions are between the lines 
cii and For small q <C 1 there are almost only stable regions except for thin resonance bands, 
while for large q > a there are almost only unstable regions, except for thin stable bands. In the 
unstable regions the characteristic exponent v will have a non-zero imaginary part. This figure is 
taken from [34]. 

3 Quantum mechanics 

Following the work of Grozdanov-Rakovic [35] and Rahav-Gilary-Fishman [31, 32] it is 
possible to use the method of separation of time scales also to calculate the quantum 
corrections to the Hamiltonian considered previously. This gives basically the effective 
Hamiltonian (2.8) (with appropriate ordering of the operators) as well as an extra quantum 
correction (to 4th order). The spirit of the calculation is to find a unitary transformation 
of the quantum Hamiltonian rendering it time independent. Again, doing so exactly is a 
highly non-trivial task but order by order in 1/co this can be done systematically and the 



-7- 



Figure 2. The frequency of the time-dependent harmonic oscillator as function of time r. 
resulting effective quantum Hamiltonian reads [31, 32] 



where x and p are now operators. Let us apply this formula to two cases in turn in the 
following. 

3.1 Quadratic potential 

In the case of Vo = A 3 x 2 the effective quantum Hamiltonian reads 

H = f + —,x 2 + 4V + 0(u^) . (3.2) 

Note that the quantum Hamiltonian to 4th order in the 1/uj expansion is exactly equal 
to its classical counterpart (eq. (2.11)) and hence also the quantum effective potential is 
confining. As before this is a special case and we can say a lot more than what is given by 
the 1/cj expansion. In fact the problem can be solved also in quantum mechanics. Lewis 
and Riesenfeld (LR) [36] found a method of invariants to solve time-dependent harmonic 
oscillator problems. They considered the generic Hamiltonian 

H = ±(p* + tf(T)x 2 ) , (3.3) 

with Q 2 (t) G R (not positive definite) and the corresponding classical equation of motion 

x + n 2 (r)x = 0, x=p. (3.4) 

Let us make a lightening review of the results of the LR-method. First we split time into 
three regions: a past (r < n) and a future (r > T2) where the frequency of the harmonic 
oscillator is constant; in between we turn on the time-dependent frequency f2(r), see fig. 2. 
If one desires, it is possible to send —t±,T2 to infinity. This problem admits the construction 
of an invariant operator called the Lewis invariant [37] 



2 



—5 + (pp — mpx) 
P 



(3.5) 



which obeys Heisenberg's equation of motion while 

1 

P 



m 2 p + n 2 (T)p- -3 =0, (3.6) 
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is the Ermakov-Pinney equation. The eigenvalues of I can be shown to be time independent 
[36], while its eigenstates are time dependent. The trick is that the eigenstates of the 
invariant operator / can be mapped to those of the Hamiltonian and hence the transition 
probabilities for any given system can be calculated in terms of the eigenstates of the 
original Hamiltonian (3.3). The general solution to eq. (3.6) for constant O is [37] 



p(r) = i^M -1 [cosh(5±sinh(5sin(2ftT/m + pj] , (3.7) 

where 5, ip £ K are real parameters. The parameter 5 encodes all information for calculating 
the transition probabilities from any given state at time r < t\ to any state at time t<i- 
Unfortunately, 5 is not analytically calculable, it can however be computed numerically. 

For our purposes, we do not really need to calculate the occupation of particular states 
in detail, we merely would like to know if the quantum system is bounded in the sense 
that the wave function yields a bounded expectation value of the position operator, or 
alternatively that the average number of created particles remains bounded. The solution 
to the Ermakov-Pinney equation (3.6) can be written as follows [38, 39] 



p(r) = \J clu 2 (t /m) + cv 2 {r/m) + 2bu(r j m)v(r j m) , (3-8) 

where u(t /m),u(r /m) are two linearly independent solutions to eq. (3.4) and ac — b 2 = 
l/W 2 , where W = uv — vii is the constant Wronskian of the solutions. 

The argument now goes as follows, if the classical system satisfies our demand on 
boundedness, then the two solutions n(r/m), v (r/m) are necessarily bounded and hence 
for any set of finite constants {a, b, c} implied by any finite boundary condition for the 
quantum auxiliary function p(r), its behavior will also obey the demand on boundedness. 
This establishes equality between boundedness of the classical and the quantum systems. 

3.2 Quartic potential 

In this case, Vo = A 5 x 4 and hence we have 

rr P 2 4A 10 fi 27A 10 , , 41 , 36^ 2 A 10 2 /nr 5n . , 

H = ^- + ^ 00 + — 3~T {?> {P» x } } + 3 4 s + O a;" 5 . 3.9 

This theory cannot be solved analytically and we have to rely on the assumption that the 
1/ui expansion is valid and equivalently the effective quantum potential is confining. 

4 Field theory 

4.1 Quadratic potential 

Consider the Lagrangian density 



£ = ±(d M 0) 2 + -m 2 cos(a;t)(/> 2 , (4.1) 



where (p(t, x) £ R is a real-valued scalar field in (<i+l)-dimensional flat spacetime, m 2 £ M>o 
and u) is an externally given fixed frequency. For convenience we choose the field to be 
tachyonic at time t = 0. The equation of motion is 

^+[a-2gcos(2r)]</> = 0, = 0(r, k) , (4.2) 
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where we have Fourier transformed the field in the spatial directions and rescaled time as 
r = ut/2 while the double-dot in eq. (4.2) denotes the second derivative with respect to r. 
Eq. (4.2) is the Mathieu equation on canonical form [40], while the parameters are related 
to the original Lagrangian density as follows 

4fc 2 2m 2 lA Q , 

a = , q = . (4.3) 

The external frequency u sets the scale of the problem, while the parameter q corresponds 
to the mass squared and finally the parameter a corresponds to the two-norm squared of 
the momentum vector, a is the parameter in the equation adding structure to the system 
due to the quantum field theory as compared to the quantum mechanics case. 

There are two types of impinging instabilities threatening the system. The first is the 
unboundedness of the potential which we have chosen to cure by making the potential 
oscillate rapidly (i.e. we consider q <C 1, see fig. 1). The second issue concerns resonant 
modes appearing upon this dynamical stabilization of the potential (i.e. for a taking on 
integer values resonant modes appear, see fig. 1). Let us emphasize that this instability is 
not that of the potential being negative but merely a secondary effect due to the oscillatory 
behavior of the potential providing opportunity for resonances in the system at hand. The 
critical momentum magnitude above which all momentum modes render the potential of 
the theory stable is 

ki > 4itical = ™ 2 • (4-4) 

Above this critical momentum only the resonances reside. For a continuum of momentum 
modes, some modes are bound to hit the resonances, which suggests that if we want stability 
we should compactify space. For simplicity we will consider only toroidal compactification 
of the spatial dimensions. Our spacetime manifold is therefore 1R x T d and we denote the 
compactification radii as Ri, i = 1, . . . , d. Now the momenta are given in terms of a set of 
integers {ni\ as follows 



h = £ , (4.5) 
Hi 

and hence 

d 4 
a = Y,0inh d i = ^2- ( 4 - 6 ) 

1 7 
1=1 1 

According to eq. (4.4) the potential becomes positive definite if we choose the compactifi- 
cation radii such that 

-^>m 2 , V* = l,...,d. (4.7) 

This ensures that there is no instability due to the potential in the theory while it does not 
guarantee the absence of resonances which can and do occur in general. In the remainder 
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of this section we will study a compactification scheme such as to avoid the resonant modes 
entirely. 

Avoiding the resonant momentum modes corresponds to choosing the parameters B{ 
such that the quantized momentum modes hit only stable zones of the phase diagram 1; 
this means they should avoid hitting any integer squared s 2 , s £ Z>o- The size of the 
unstable band can be estimated for large s as [34] 



A, 



(4.8) 



which also calls for very small q <C 1 • The stability condition now amounts to the following 
inequality 



i=l 



> A s , 



(4: 



which we will simplify by setting all the 9s equal 9i = 9 (hence Ri = R) and thus 



d 2 
nf 



i=l 



-sA 



> 



V>(«) 



(4.10) 



where A = — log q > 0. 

Instead of addressing this problem directly, let us consider an easier but as we shall 
see, more general problem, viz. consider the stability condition 

1 



> Hi) , 



(4.11) 



where p,q S Z>o are integers. Now the strategy is to find the best approximation p/q to 
1/9 and try to estimate how big their difference will be as function of q. If it will be big 
enough, i.e. bigger than tp(q), then we are guaranteed absence of resonances. 
Let us define 9' = 1/9 and consider the Diophantine approximation to 9' 



P 



< 



1 



Mq 2 ' 



(4.12) 



where M is a constant of order one depending on 9'. In general 9' is a real number but we 
will consider the case where it is an irrational number, 9' S R\Q and positive. It will prove 
convenient to utilize the continued fraction as a method for obtaining the Diophantine 
approximation p/q to the irrational number 9' . The continued fraction is given by 

1 



Pn 

— = a H 

Qn ai + 



(4.13) 



a 2 + 



a n -i + 



where {ai} are positive integers and p n /q n is called the n-th principal convergent, while it 
holds that 



lim 



(4.14) 
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Hence if 9' is rational a solution exists. On the other hand, if 9' is irrational, there will 
only exist an arbitrarily good approximation in form of a rational number. The question is 
how fast it converges, viz. how good the rational number pjq approximates 9' as function 
of q. We need to invoke a few theorems due to Euler and Lagrange in order to answer that 
question. 

Theorem 1 [41, 42] For n > we have 

q n p n -l -PnQn-l = (-1)™ ■ (4.15) 

Corollary 1 [4-1, 42] If a±,a2, ■ ■ ■ are positive integers, then p n and q n are relatively prime 
and < qi < q% < ■ ■ ■ forms a strictly increasing sequence of integers. 

Theorem 2 [41, 42] The best approximations to 9' are the principal convergents to 9' . In 
fact for n > 1, q n is the smallest integer q > q n -\ such that \ \q9'\\ < \\q n ^i9'\\. 2 

In view of the above theorem and the fact that {q n } is a strictly increasing sequence of 
integers, we just need to determine how good the continued fraction approximates 9' for a 
given q n . 

Theorem 3 [41, 42] For even n, the n-th principal convergents of 9' form a strictly in- 
creasing sequence converging to 9'. For odd n, the n-th principal convergents of 9' form a 
strictly decreasing sequence converging to 9' . Furthermore, we have 

-—< <\q n 0'-Pn\< • (4.16) 

^Qn+l Qn+1 + Qn Qn+1 



It then follows from theorem 3 that 



]' _ En 

Qn 



> 



(4.17) 



which gives us an estimate on how bad the approximation is iff we know how much larger 
q n +i will be compared to q n . It can be shown that 



Qn+l — an+lQn + Qn-1 , 

holds. Using corollary 1 we can infer that 

Qn+l < ifln+l + 1) 

and hence we can write eq. (4.17) as 



Ih, 



> 



1 



2(a n+1 + l)g2 



(4.18) 



(4.19) 



(4.20) 



We adopt the notation ||£|| as the distance between £ and the nearest integer. 
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Comparing the two right-hand sides of eqs. (4.20) and (4.11) we arrive at the following 
condition 

Q 

a n+1 + 1 < -q^- 1 exp{q n A) , (4.21) 

which when satisfied guaranties absence of resonances. 

An example which satisfies the condition (4.21) is the golden ratio: 1/9 = which 
has the integers Oj = l,i = 1,2,... and is therefore the slowest converging irrational 
number. Using eq. (4.18) it is seen the q n s form a Fibonacci series and hence q n < q n +i < 
2q n+ i. Plugging the a^s into eq. (4.21) the left-hand side is constant while the right-hand 
side is rapidly growing with re. Another example is 1/9 = e having cyclic a^s: [ao, ax, ■ ■ ■ ] = 
[2, 1, 2, 1, 1, 4, 1, 1, 6, 1, 1, 8, 1, 1, 10, • • • ]. This number also satisfies the inequality (4.21). 
Square roots of integers always have a periodic series in {a^} and hence are easy to check 
if they satisfy eq. (4.21). 

The condition (4.21) is valid for eq. (4.11) which is a generalization of the problem we 
want to solve, namely that of eq. (4.10). Yli=i n i/ S ' 2 * s n °t an arbitrary rational number 
p/q but (for finite dimensions d < oo) it is a subset of rational numbers, so the condition 
(4.21) is stronger than what we need. Thus if it is satisfied also eq. (4.10) is satisfied and 
hence there will not exist any momentum mode in the compactified field theory hitting the 
instability bands of the phase diagram of the Mathieu equation (no resonant modes). 

4.2 Damage control 

Let us consider how bad things go if one should not take our advise on the suitable values 
of 9. We will here consider the consequences in the case that the chosen value of 9 gives 
rise to one or more modes hitting a resonance band. We can estimate the time t\ it takes 
until the norm of the unstable mode has reached the size A as 

2 R\/9 
w|<y(i/)| I^H 17 )! 

In addition we need an estimate of the size of the imaginary part of the characteristic 
exponent v, which can be found in the literature [43], valid for q a (an alternative 
method for estimating v can be found in app. A) 



~ / x \/ (o-s — a)(a — b s ) 

9(i/)~±-¥-^ — ^ -, (4.23) 

where a is the parameter of the Mathieu equation, while a s (b s ) is the upper (lower) 
boundary of the s-th instability band. Let us assume the worst case scenario, i.e. the mode 
hitting right in the middle of the instability band giving the maximal imaginary part of v 
in which case (for g<l) o= (a s + b s )/2 and hence 

~ £ = £ , (4-24) 

where A s is the size of the s-th instability band (4.8). This value is a rapidly decreasing 
function of s and hence the "lifetime" in eq. (4.22) will be parametrically large for high 
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enough s. Now that one can tolerate the large life times associated with large values of s, 
one needs to remove the thread of resonances for smaller values of s. There are two tools 
available for that; an appropriate choice of 6 or a sufficiently small Ruj. 

4.3 Quantizing the compactified field theory 

We have now obtained a classical field theory which under certain conditions is stable 
in the sense that for each momentum mode, the equation of motion gives a - from the 
other modes decoupled - bounded wave function under time evolution. For a single time 
dependent harmonic oscillator we have seen in sec. 3.1 that the quantum state evolves in 
a controlled manner basically determined by the classical behavior. Let us remark that 
Q 2 (r) of eq. (3.4) is allowed to contain a non-zero constant term a which is appearing in 
eq. (4.2). The time evolution of the quantum field theory state can be explicitly calculated 
[44] which is basically done by acting with the Hamiltonian on a state vector living in a 
Schwartz space in the Schrodinger representation and finally composing an infinite product 
of these over all momentum modes. With respect to our stability/boundedness criteria, we 
are basically home free due to the fact that we are dealing with a free (i.e. non-interacting) 
field theory whose infinitely many harmonic oscillators are mutually uncoupled. We would 
however like to address the question of whether the average number of particles created 
remains finite, summing over all momentum modes. For that we use the solution (3.8) 
subject to the boundary conditions p(0) = {On 2 — 2q)~ 1 / 4 and p{0) = which in turn we 
choose to feed with the Mathieu cosine and sine functions expanded in small q (see app. B). 
For the characteristic exponent v we use the approximation (B.6). In order to calculate 
the transition probability we need a quantum mechanical integration constant 5 calculated 
for each mode n. The following integrated form of eq. (3.6) [36] (for m = 1) 



p 2 + fi 2 p 2 + -o =2101 cosh <5 
P 2 



n 2 = 6n 2 - 2gcos(2r) 



(4.25) 



defines 5. Plugging in the above described solution p we find 





- 2(1 + \/0n)cos 2(1 -\f6n)T - 2(1 - v^n) cos 2(1 + \/!)n)T } + 0(q 3 ) , 



with r > 



The average number of particles N [36] is then 




(4.27) 



n 1: ...,n d =l 




Expanding the above expression in large n>lwe obtain the 



average number of created particles 




oo 



(4.28) 
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from which it is seen that (N) converges for d < 5 spatial dimensions. Similarly we can 
calculate the total energy of the system, neglecting the mass of the created particles, as 

Eiot -R E I (cosh (5-1), (4.29) 

m,...,n d =l 

which then converges for d < 4 spatial dimensions. 

4.4 Quartic potential 

For the field theory 

^ = ^(^) 2 -^cosM)0 4 , (4.30) 

the detailed information which was available for the quadratic potential is not readily acces- 
sible. In the following we will take a leap of faith and assume that the system thermalizes 
due to the fact that it is an interacting theory. The actual calculations we will perform 
will be related to the ABJM theory [45] on a sphere. In the next section we consider the 
gravitational dual of such a theory with time-dependent couplings. In the process we will 
provide evidence from the bulk that the theory indeed thermalizes. 



5 AdS/CFT 

We will study the gravitational bulk manifestations in the regions of (in-)stabilities obtained 
for the boundary field theory. Our expectations are that the gravitational dual of an 
unstable boundary potential will result in a big crunch. For those values of parameters for 
which the boundary theory is eventually stable, we expect that a black hole (BH) horizon 
will form. We will consider the setup of 1 1-dimensional supergravity in AdS4 x (S /%k) as 
a ground on which to identify the appropriate bulk field configurations, corresponding to 
the different types of boundary potentials we study. N = 8 supergravity in four dimensions 
describes the massless sector of the compactified 1 1-dimensional supergravity on S 7 . It is 
possible to introduce a consistent truncation involving only gravity as well as a single scalar 
field 4> [46]. We thus consider the following effective action 



(|-^OT-y(0) 



with V(4>) = -2 - cosh(\/20) « -3 - 4> 2 + 0(cp 4 ). As pointed out in [7, 8], the same 
consistent truncation can be used for the more general compactification of 1 1-dimensional 
supergravity on the quotient space S /Zfc. The scalar field (j) has a mass squared m 2 = —2, 
which is above the Breitenlohner-Freedman bound [47] wig F = —9/4. The AdS metric in 
global coordinates can be written as 



ds z = _ (1 + r 2) dt 2 + Jj_ + r 2 dQ 2 = _Jt + + tan 2 x d tf (5 .2) 

1 + r z cos z x cos z x 
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where r = tan cc. The scalar field goes to zero asymptotically as 

^)«" + ^- ( 5 - 3 ) 

According to the AdS/CFT correspondence, M-theory in global AdS4 x (S 7 /Zj.) is dual 
to ABJM theory [45] on S 2 x R, i.e. a U(A) x XJ(N) superconformal Chern-Simons (CS) 
theory with matter and CS levels +k and —k. It contains four M = 2 chiral superfields 
Y s transforming in the (JV, iV) representation of the gauge group. Due to the conformal 
coupling, the four scalars Y s acquire masses proportional to I/-R52 which we set to unity. 
The field <fi of the consistent truncation in eq. (5.1) corresponds to the following scalar mass 
operator in ABJM theory 

O = ^Tr (YiY} + Y 2 Y\ - Y z Yl - Y 4 Y^ . (5.4) 

Multitrace deformations [26, 27] of the form S = So — N 2 W(0), correspond to the following 
boundary condition 

fi = W'(a) . (5.5) 
This is possible due to the fact that the mass of the scalar field lies in the interval 

- - < m z < -- , (5.6) 
4 4 

and hence both scalar-wave fall offs of eq. (5.3) are normalizable. 

By appropriately choosing W(0), it is possible to construct unbounded potentials in 
the boundary theory for which Y s = is a metastable vacuum; for instance by choosing 
W = H^0 3 /3 with an arbitrary sign for H3, or alternatively W = H2O 2 /2 with H2 < 0. 
The first case corresponds to AdS-invariant boundary conditions [48]. We first study the 
case of a time-independent unstable potential and identify the bulk signature of that and 
then we introduce time dependence [2]. 

5.1 Instantons and crunches 

Let us first consider the case of a time-independent multitrace operator. When the multi- 
trace operator renders the boundary theory metastable, this metastability manifests itself 
in the bulk as a Coleman-de Luccia instanton [1] when the multitrace operator is marginal 
and by some other type of bubble when it is a relevant operator [10, 13]. In the AdS4 setting 
we are considering, instantons were studied in detail in [2, 3]. The following spherically 
symmetric Ansatz for the metric is used 

ds2 = {m +p2dnl )- (5 - 7) 

The Euclidean action reads 

S = Jd 4 x^ f-| + l -d^4> + V^U - J 9 d z x^h c K + S ct + Sg , (5.8) 
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where h c is the induced metric on the cutoff surface, K is the trace of the extrinsic curvature 
of the boundary (giving the Gibbons-Hawking term), S ct and Sq contain counter terms and 
the boundary terms, respectively. In the AdS4 setup with m? = —2 which we study, S, 
and Sq take the forms [49, 50] 



ct 



S 9 = f d 3 x Vh (W(a) - aW'{a) + a/5) , 
Jd 



where A_ = 1, h is the induced metric on the boundary and R{h c ) is the Ricci scalar of 
the induced metric on the cutoff surface. 

Using the p coordinate, the equation of motion for the instanton is 



P 



(5.9) 



For numerical convenience we make a change of variables as follows; instead of p we will use 
x = tan -1 p; and instead of b, A = b 2 /{p 2 + 1). The equation of motion for the instanton 
then reads 



+ 



2 + cos 2 x A' 
-. + 



sin x cos x 2 A J A cos 2 x 
which is sometimes conveniently written as 



0. 



,4 



2V{4>) sin 2 x — 6 cos 2 x 
sin 2 x cos 2 x{(f)') 2 — 6 



(5.10) 



1 



V sin2x „ r . ,. 



COS' X 



6 



2 + cos 2 x . sin 2x 



12 



+ 



sm x cos x 



(2 + sin 2 x) {4)'f . (5.11) 



Fig. 3 shows an example of an instanton solution. The family of solutions is parametrized 
by the value 0o — ^(0), i.e. the value of the scalar field at the origin. Each of these 
solutions corresponds to a particular value of (a, 0) of eq. (5.3) at infinity; see fig. 4. The 
limit (f)Q — > corresponds to #3 —> ±00 for tripletrace deformations, while for doubletrace 
deformations H2 — > — 1. 





(a) <j>(a 



(b) A(a 



Figure 3. Instanton solution with <f)Q = 0.1. 
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Figure 4. Left panel: a (solid line) and (3 (dashed line) for the instanton as function of 4>q. Right 
panel: i/ 2 = fi/ot as function of <f)Q. 

The decay rate can be computed in terms of the Euclidean action of the instanton, 
using the Coleman-de Luccia [1] method. The relevant quantity determining the tunneling 
rate is e~( s ~ SAdS ^ h , with S being the instanton action and S^s the- void AdS action. In 
the tripletrace case, it was checked in [2] that 5 — SAdS i s a finite quantity and hence the 
tunneling occurs in a finite time. We checked that this is also the case for the doubletrace 
deformation with negative H2. 

After the tunneling occurs, the initial condition is obtained by restricting the instanton 
solution to the equator. In the case of AdS-invariant boundary conditions with a tripletrace 
deformation, time evolution can be obtained directly via analytical continuation of the 
instanton. Inside the lightcone the metric has the form of a Friedmann-Robertson- Walker 
solution 

ds 2 = -dt 2 + a 2 (t) (dr 2 + sinh 2 r dttf) , (5.12) 
with equations of motion 

« 2 = 1 + y ( y + v(m , 4> + £j> + v'W = , (5.13) 

describing a universe which is initially expanding but then eventually contracting; thus 
developing a big crunch singularity within a finite time [1]. 

5.2 Black holes 

We recall that the Einstein-scalar theory (5.1) has also other time-independent solutions 
than empty AdS. For all values of H2 and H3 the AdS-Schwarzschild BH is a solution with 
the metric 



Ae 25 , 2 dx 2 o ,^2 a ^bh cos 2 x 

=—dt 2 + - t— + tan 2 x dm , A = l 

cos^ x A cos z x 47r tan x 



with (5 = 0. The mass and the temperature as a function of the horizon Xh are 

tanrr/j If 1 \ 

M B h = 47r — ^ — , T B h = ~r 3 tan x h + , (5.15) 

cos z Xh 47r \ tan x^ 
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and the specific heat is negative for Xh < tt/6. 

More general BH solutions with scalar hair have been studied in [48]; among the 
time-independent solutions, there are also solitons without horizons, see e.g. [2, 51]. The 
time-independent equations of motion read 

1 / 9 a -fi-r\' dV x 1 „ r , sin x cos x T? .„ ^ . 

tan 2 xAe d ^ - — e~ b =— = 0, 6' = * 2 , (5.16) 



tan 2 x V / dip cos 2 x 2 

l + 2sin 2 x. ,. sinxcosx^n , rr , . /- -,-\ 

A = 1 - A) Al' 2 - + 3 tana; , 5.17 

sin x cos x 2 

where ^ = d<j)/dx. At the horizon x^ (defined by A{xj l ) = 0) the following boundary 
condition should be imposed 

tanx 1 dV , 

^ = . (5.18) 

1 — V((f>) tan 2 x cos 2 x dcf) 

All these time-independent solutions [2, 48, 51] have a negative ratio H2 = (3/a < 0. The 
only time-independent solution compatible with H2 > is the Schwarzschild BH, which is 
a solution for arbitrary Hi- 

We will next proceed to uncovering the duals of time-dependent boundary theories; 
the BH will be a good candidate to screen the big-crunch singularity for those cases where 
we will find that dynamical stabilization occurs. 

5.3 Time-dependent equations of motion 

We now turn to matching the gravitational solutions to time-dependent boundary condi- 
tions. This requires in turn, that the bulk solutions be time dependent. To this end we 
consider a spherically symmetric setup, with the metric 

e -M(x,t)A(x,t) l2 dx 2 2 

ds 2 = V ' ! dt 2 + — r ^ + tan 2 xc&7 2 .. (5.19) 

cos^ x A[x, t) cos z x 

as well as a scalar field <p(x,t). The cross term dxdt is absent from the metric due to the 
choice of gauge. There is still some freedom in shifting the variable 5 by a function of time; 
this corresponds a time re-parametrization. In the following we will set 5 = at x = tt/2, 
such that t is the global time coordinate in asymptotic AdS space. 

Let us introduce the auxiliary variables VP = (f>' and n = A~ 1 e s <j) (where ' denotes 
derivative with respect x and ' with respect to t) for which the equations of motion read 

(Ae- s n)' = V, -4r- (tan 2 xAe-^V - -^e~ s —^- = tl, (5.20) 
tarn x V / d(p cos z x 

l + 2sin 2 x. sin x cos x ,. t9 :S 

A = (1 -A) AC& + n 2 - (V(6) + 3 tan x , 5.21 

sm x cos x 2 

„/ Sill X COS X ,^^0 T o, ,/ _ , v 

5' = (n 2 + ^ 2 ) , = (5.22) 

<p = UAe~ 5 , A + (sinxcosx)A 2 e _<5 ^n = 0. (5.23) 

The last equation of (5.23) is a consequence of the other ones. 
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In order for the solution to be smooth, we must require that \J/| x= o = n'| x= o = ^4'| x =o — 
n U^7r/2 = A '\x^n/2 = and A\ x=0 = A\ x ^ n/2 = 1; both n| x=0 and ^\ x ^ n/2 are non- 
vanishing functions of time. The quantities (a, f3) are directly related to the values of \P 
and if' at the boundary 

(tt/2) = -a , *'(tt/ 2 ) = 2/3 . (5.24) 

Near the boundary, </> is going to zero linearly in (tt/2 — x), while is finite. In the following 
we will specialize our treatment to time-dependent boundary conditions of the form 

W = ^0\ H(t) = H 2 (t) = ^ (5.25) 
2 a{t) 

corresponding to doubletrace operators. The boundary conditions at x — > ir/2 are 

= , IT 7 = -2HV - 2HW . (5.26) 

5.4 Observables 

In this section we will define some useful observables describing geometric features of the 
numerical solutions. The first observable we address is the scalar curvature, which by 
means of the equations of motion, is 

R = W(<p) +^(^ 2 -n 2 ) cos 2 x. (5.27) 

The second observable we need is the total energy for which it proves convenient to 
express the metric component g rr in terms of the function A(x). Setting A(ir/2) = 1 and 
A! '(vr/2) = we get 

1 / A"{ir/2) \ 1 Ajgfr/g) 1 
^ = ^-^ 1 + — 2~ + 6 ^ + '-- (5 ' 28) 

Expanding eq. (5.21) in series around x = 7r/2, it follows that A"(ir/2) = ^ 2 (ir/2) = a 2 . 
Using an expression for the total energy given in [51], we obtain 

U^(7r/2) , ^ , \ A _(AV>(*/2) , 3 ^ 



M = Att y + af3 + J P(a)daj = 4tt y + -i?a ^ J . (5.29) 

This quantity is in general not constant, because the time-dependent multitrace operator 
acts as an external source of energy for the conformal field theory described by the AdS4 
geometry. 

A trapped two-surface S (see e.g. [52]) is defined as a closed surface with the property 
that the expansion scalars 

Oi = (g ab + l a n b + l b n a )VJ b , n = (g ab + l a n b + l b n a )V a n b . (5.30) 

in each of the two forward-in-time null directions Z a) ra b (normal to S) are both negative. 
We use the cross normalization l a n a = —1. Then, considering a foliation by spacelike 
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three-surfaces a point q G £j is said to be trapped if it lies on a trapped two-surface 
5 in T,f. The apparent horizon in £4 is thus constructed as the boundary of the union of 
all the trapped points. Let us consider the following null direction of ingoing and outgoing 
wave fronts in the (t, x) direction 

l a = -—l {- A ,-e s ,0,0}, n a = -=± {-A,e 5 ,0,0}. (5.31) 

\/2Ae° cos x \/2Ae° cos x 

An explicit calculation gives 

\f2A 

0l,n = T-. • (5.32) 

smx 

This shows that no apparent horizon forms as long as A remains positive. 

The metric (5.19) is an example of polar coordinates, which are in general defined by 
the condition Kg = K,t = 0, where K, is the extrinsic curvature tensor. It is a general 
feature of polar coordinates that they can not be used to parametrize events happening 
inside apparent horizons (see e.g. [52, 53] for a discussion); instead they become singular 
when an apparent horizon forms. Apparent horizons are then detected by the condition 
that the metric (5.19) is singular [52, 53], which is therefore given by a vanishing A. An 
alternative foliation and metric Ansatz should be used in order to penetrate the apparent 
horizons. 

5.5 Numerical results 

In this section we present numerical solutions of eqs. (5.20-5.23) obtained for different time- 
dependent boundary conditions, specified by the function H(t) in eq. (5.25). As initial 
condition we use the instanton obtained by eqs. (5.10,5.11) restricted to the equator. We 
imagine that H (t) is constant for t < and that the system is in a metastable vacuum. We 
then assume that a tunneling occurs at t = from where we introduce time dependence to 
the doubletrace deformations. As a parametrization of the instanton solutions, we fix (fro, 
corresponding to the value of the scalar field at the origin of the instanton. The relation 
between </>o and the initial value of H(t) = H2 is shown in fig. 4. 

We first consider the case of constant H(t) = ho < for which the solution corresponds 
to big crunch geometries. A curvature singularity (which extends to the boundary) develops 
in a finite time. The energy can be computed via eq. (5.29) which we have checked is 
diverging at the time when the crunch singularity forms. A plot of the time it takes for 
each instanton to evolve into a crunch as function of 4>o is shown in fig. 5. This quantity 
increases with the steepness of the potential. The reason for this is that the larger the 
value of (fio is, the larger is the energy which causes the crunch. Note on the other hand 
the smaller the value of <fio is, the faster is the decay of the vacuum. This behavior can be 
seen in the right-hand part of fig. 4. 

Let us next consider a boundary condition of the type 

H(t) = h (2e~ t2 - l) , ho < , (5.33) 

which from the boundary point of view corresponds to a negative and unbounded potential 
at the initial time, but is then lifted in a finite time ~ 2 to a positive definite one. A 
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Figure 5. The time it takes for an instanton to evolve into a big crunch singularity for constant 
and negative H(t) = ho < as function of 4>q. 



solution within this class is shown in fig. 6. For a massless particle in AdS space, it takes 
the time n to reach the boundary and come back to the origin, corresponding to a frequency 
wo = 2. Indeed in the plots both the fields and the curvature oscillate with a frequency of 
roughly 2. 

At late times we expect a BH to form due to the weak turbulent instability discussed 
in [22, 54, 55]. However, in the regime of small (j>o the horizon forms only after a very long 
time and after the scalar field has bounced off the AdS boundary many times. The time 
of BH formation as function of 4>q is shown in fig. 7. For computational constraints we 
computed the solutions only till time i max — 40. In the cases where a horizon has formed 
during the finite calculation time, we have checked that the final configuration is, to a good 
approximation, a Schwarzschild BH outside the horizon. 

We are now ready to study the boundary conditions related to dynamical stabilization. 
These alternate between a negative, unbounded potential and a positive one 

H(t) = h cos(ujt) , h < . (5.34) 

For sufficiently small u a crunch will occur, otherwise given a general point in the parameter 
space ((j>o,tj) a crunch may or may not form at some point in the future. A plot of the 
parameters for which a big crunch singularity forms (before time i max = 40) is shown in 
fig. 8. 

The transition between the crunch and the BH regime should correspond to a BH with 
infinite mass. This regime however is challenging to study numerically as it is sensitive to 
the UV cutoff 

%m&x which in turn requires increasingly small finite difference steps. 
An example of a solution where a horizon forms is shown in fig. 9. As explained in 
section 5.4, a horizon can be detected by virtue of a vanishing A. The scalar curvature 
R becomes quite large inside the horizon. The time evolution after the time where A 




(c) 6(x,t) (d) scalar curvature R(x,t) 



Figure 6. A solution with boundary condition (5.33) and (j>o = 0.1. The initial negative unbounded 
potential is lifted to a constant positive definite one in a finite time t ~ 2 much before the time when 
the crunch singularity would have formed if the unbounded potential were kept constant (see figure 
5). The resulting configuration in this case is just a scalar wave reflected by the AdS boundary 
conditions. 

vanishes cannot be trusted, because the metric becomes singular and hence cannot be used 
to parametrize events inside the apparent horizon. 

An example of a parameter point where no horizon forms before time i max = 40 is 
shown in fig. 10. Here the scalar curvature remains rather small everywhere. The curvature 
and the scalar field at x = are shown as function of time in fig. 11. 

For sufficiently large frequencies, we have checked that in the case that no horizon has 
formed before i max , the maximum of the energy and the maximum of \R\ remained small 
(see figure 12). When a horizon forms, the scalar curvature tends to become large; 
however this all happens inside of the horizon. The energy of the resulting BH does not 
diverge, on the contrary it tends to become smaller for increasingly larger frequencies. In 
fig. 13 is shown the time of horizon formation as function of ui for (f>Q = 0.1. 

There are sharp transitions in the time of BH formation as a function of u (see fig. 13). 
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Figure 7. The time after which a BH is formed as function of </>o, for the boundary condition 
(5.33). The simulation is run until t max = 40. For </>o < 0.46 we notice a very sharp transition: no 
horizon forms before t max ; for <fio = 0.46, it is found that a BH horizon forms at t ps 882. 
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Figure 8. Shown in light gray, points in parameter space giving rise to a crunch before i max = 40, 
for the boundary condition in eq. (5.34). The points for which a BH is formed before t max are shown 
in dark gray. Points in parameter space for which neither a crunch nor a BH is formed before i max 
are shown in white. 

For frequencies u around even integer values, a BH forms rather quickly, whereas for 
frequencies u near odd integer values, we have observed no BH formation before stopping 
the numerical code at t = 800. A field propagating with the speed of light takes the time 
7r to reach the boundary and return from x = 0; this corresponds exactly to the frequency 
ujq = 2. Perhaps this provides some qualitative explanation for how come BHs form more 
easily subject to even frequencies. 

Two options are possible upon horizon formation. Either the solution tends to a 
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(a) <j){x,t) 



(b) A(x,t) 




Figure 9. A solution where a BH forms for 0o = 0.1 and us = 0.9, corresponding to a dark gray 
point in fig. 8. Dynamical stabilization is at work: the big crunch singularity does not reach the 
boundary in a finite time as it is screened by an apparent horizon, which is where A = 0. The 
absolute value of the scalar curvature (d) becomes large (signaling the formation of a singularity), 
but it happens inside the horizon. 

Schwarzschild BH at late times, as it is the only time-independent solution compatible 
with our boundary conditions or alternatively the BH will be time-dependent and its size 
and mass will increase indefinitely - we are effectively pumping energy into the system. 
Unfortunately our choice of metric Ansatz (5.19) is not permitting us to study this issue, 
because it becomes singular as soon as the horizon has formed. 

In the case of a massless scalar in flat 4-dimensional spacetime, Choptuik found for 
an arbitrary family of solutions parametrized by p, that both the critical horizon radius 
and the BH mass scale as (p — p*) 7 , where 7 0.37 and p* is the critical value for BH 
formation [56]. An identical scaling was found for AdS4 with a massless scalar in [22, 57]. 
For a massive field in flat space (e.g. [58]), depending on the class of solutions considered, 
it is also possible that BH formation begins at finite mass and radius (type I transition), 
instead of beginning at infinitesimal mass (type II transition). In the case of a type II 




(c) 5(x,t) (d) scalar curvature R(x,t) 



Figure 10. A solution where no BH forms before i max = 40 for fa = 0.1 and uj = 3.0, corre- 
sponding to a white point in fig. 8. Dynamical stabilization is also at work here, even if no BH 
horizon has formed before i ma x- The solution corresponds to a scalar wave being reflected many 
times by the time-dependent boundary condition, without a dramatic increase of the total energy. 
The field configuration is similar to that of fig. 6, even though the stabilization in the latter case is 
not dynamical. 

transition, the same critical exponent 7 was found as in the massless case. The exponent 
7 is not universal for all kinds of matter; for example, in the case of type II transitions in 
Einstein- Yang-Mills equations [59], 7 ~ 0.2, (for a review of critical gravitational collapse 
with various kinds of matter, see [60]). 

In AdS4 with a tachyonic scalar in the mass range (5.6), it is in principle possible to 
study both the critical behavior as function of the initial field profile, as well as function of 
the boundary condition at x — > tt/2. The boundary condition (5.33) gives rise to a type II 
transition, which we have observed due to the formation of arbitrarily small BH horizons 
with critical value (fio c ~ 0.468. On the other hand, for the boundary conditions (5.34) 
and fixed <pQ = 0.1, there are several type I transitions as a function of co: for example, 
BH formation begins at radius ~ 0.5 and t = 64 for oj = 1.175 > u) c , while no BH horizon 
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Figure 11. The value of 4> and the scalar curvature R at x = are shown as functions of time for 
u> = 3 and 0o = 0.1. For comparison also the driving potential H (t) is shown. 




Figure 12. Left panel: Maximum of the energy as function of ui for <f>^ = 0.1. Right panel: 
Maximum and minimum of the scalar curvature R as functions of w. The *s correspond to the 
minimum of A in some units, showing the correlation between the latter quantities and the formation 
of a horizon. 

has formed before t = 800 for lj = 1.17 < uj c . We leave a more accurate study of these 
transitions and a determination of the critical exponents in case of a type II transition for 
future work. 

5.6 Numerical scheme 

We solve the time-dependent PDEs using the leapfrog method with a five-point - 4th order 
difference stencil for the spatial derivatives, a grid of 2500 points and Courant number 0.5. 
Explicit dissipative terms are added. The instanton solution (taken as the initial condition) 
is obtained by solving eqs. (5.10,5.11) using 5th order Runge-Kutta. In order to initialize 
the leapfrog algorithm, we need to provide the first and second order time derivatives of 
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Figure 13. The time of horizon formation as function of u for <po = 0.1. In some intervals of 
frequency centered around w taking on even integer values, a BH horizon formation occurs rather 
quickly. On the other hand, we have observed no BH formation for odd integer values, u = 3, 5, 7, 
in calculations running up to the time t = 800. Horizon formation at t = 220 is observed for w = 1, 
while no BH formation is observed up to t — 800 for ui = 1.17. 

the initial instanton profile function. We set all the first order time derivatives to zero, 
except 

n = — . (5.35) 

tanx 

The second derivatives read 

\ cos x sin x y tan z x J J tan x 

1 = A 2 e~ 25 ^ 2 cos 2 x , (5.36) 

and II = 0. No simple expression for <5 is available and hence 5 is recalculated using 
eq. (5.22). For practical reasons we introduce an IR and a UV cutoff x m \ n = 0.001 and 
Xmax = 1-565 (if we use values of x max closer to 7r/2, the corresponding boundary conditions 
make the numerics unstable). In order to update c/> at the boundary x max and A at the 
origin x m \ a we use eqs. (5.23). Subsequently eqs. (5.21,5.22) are solved using a 5th order 
Runge-Kutta method. The boundary value <5(x max ) is fixed to be zero. The eqs. (5.23) are 
then used to estimate the numerical accuracy. 

6 Discussion 

In this note we have discussed the issue of dynamical stability in various examples of theories 
with time-dependent potentials which are periodically unbounded from below. Considering 
an oscillating quadratic potential in the case of classical and quantum mechanics as well 
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as in field theory, the stability regions are encoded by the phase diagram of the Mathieu 
equation. The zero mode is stabilized when the frequency is high enough. In addition we 
found a compactification scheme relying on number theory such as to avoid resonances in 
all momentum modes of the field theory. The effective Hamiltonian for periodically driven 
systems [31, 32] suggests that a similar dynamical stabilization occurs also for interacting 
theories. We have found numerical evidence for this in a theory with an AdS/CFT dual; for 
sufficiently large u the big-crunch singularity is screened by a black hole horizon, preventing 
the crunch singularity to reach the boundary in a finite time. We have observed sharp 
transitions in the black hole formation time as function of the frequency uj. For frequencies 
near even integer values, horizons form rather quickly while for those near odd integer 
values, we have observed black hole formation at very late times or not at all, due to finite 
simulation time. 

In fact we would be curious to see if such a potential could be a holographic description 
of a universe oscillating between bangs and crunches or healed bangs and healed crunches. 
An analytic study of such system was carried out in [61-65] and references therein, but there 
the system had an infinite amount of boundaries and in the presence of closed time-like 
curves. One could wish that all this could be encapsulated on a single boundary without 
redundant features. This could have been a system for which the time-dependent couplings 
on the boundary are appropriate. Unfortunately, for this, the numerical analysis is not a 
substitute for an exact solution. It is worthwhile to study such systems non-perturbatively 
even if there is a lack of a strong phenomenological motivation at this stage. 

Let us comment on possible topics of future directions: 

• For the free field theory part one could consider different compactifications of the 
spatial dimensions to see if this gives rise to additional structure. 

• It would be challenging to follow the black hole evolution after the horizon has formed; 
for this it might be necessary to use a different (non-polar) metric Ansatz and also 
introduce a singularity excision, e.g. as in [21, 66, 67]. 

• We did not address here the precise nature of the transition between the big crunch 
and black hole of infinite mass; for example, it would be interesting to determine the 
dependence of the black hole mass near the critical frequency oj c separating the black 
hole regime from the crunch. 

• A more detailed study of the transition type (I or II) at the threshold of black hole 
formation, both as a function of initial field profile (fio and of the time-dependent 
boundary condition uj, would be worthwhile to investigate; in particular calculating 
the critical exponents for the type II transition. 
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A Estimating the characteristic exponent 

The characteristic exponent v can be estimated numerically by calculating the solution 
y(r) to eq. (4.2) subject to the boundary conditions y(0) = 1, y(0) = 0, and then evaluated 
at r = vr [34] 

cos(^7r) = y(ir) , (A.l) 
where v is implicitly a function of a, q. 



B Estimating the amplitude 

From the interpretation that the generalized coordinate x describes a particle, a natural 
question would be how far the particle goes. Let us first consider the simple example of a 
particle at rest and situated at position x(0) = 1 at time r = 0. Using an expansion valid 
for small g<Cl and v real but not an integer, the solution reads [34] 



r I qcos[(z/+2)r] qcos[(^-2)r] q 2 cos[Q+4)-r] q 2 cos[(^-4)r] q 2 (v 2 + l) cosQt) 
COS[VT\ I ^ TT ^ 1 h 32(^+1)^+2) + 32(^-1)^-2) 16(^-1)^ 



1 _ S _L 2 J 2! 1 2! q 2 (y 2 +l) 

1 A{u+1) 4(i/-l) ~r 32(y+l)(v+2) 32(i/-l)(i/-2) l&(v 2 -l) 2 

0(q 3 ), (B.l) 



which is called the Mathieu cosine function and from this we can estimate the amplitude 
max(|C|) = 1 + i + - f** „ + C% 3 ) , (B.2) 



a = 1 - v 2 , /3 = 5 + v 2 , for 1/ < 1 , 

a = 2 + 2i/ , £ = - 11+ l 5 { ^2) 2+u3 , for K ^ < 2 , 
a = 2 + 2zy , /3 = 3 - 2v + v 2 , for v > 2 . 

Let us first consider a = 0. Using 

. , i- \ irq 2 siniJair) 

cos(i/7r) ~ cos(VaTr) + : f=7 tt - , (B.3) 

4va(a - 1) 



we obtain 



v-J^ + Otf), (B.4) 



which upon insertion in the amplitude (B.2) gives 

Bq 



2 



max(|C|) = l + q+-^ + 0(q d ). (B.5) 
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For a > 0, the expansion of v changes according to 

v~y/Z-(-l)WSl f (B.6) 
4^a(a - 1) 

which to order q 2 then gives rise to the amplitude (B.2) with v = yfa. The above amplitudes 
tell us how far the particle will go, i.e. the absolute value of the largest value that the 
coordinate \x{€)\ will take on, with the initial condition that the particle starts at rest at 
x(0) = 1 and time r = 0, assuming that (a, q) are in a stability band of the phase diagram, 
see fig. 1. If one wishes to start with the particle situated at x$, it suffices to multiply the 
solution (B.l) by xq since the equation of motion is linear and hence also the amplitude 
scales accordingly. 

A similar situation can be estimated as well, namely considering the particle situated 
at the origin x(0) = but having velocity x(0) = 1 at time r = 



S(r) 



oItJt/tI - gsin[(^+2)r] qsin[(t/-2)r] q 2 sin[(i;+4)T] q 2 sin[(^-4)r] _ q 2 {u 2 + l) smjur) 

SlIl^Tj 4 (H-1) + 4(i/-l) 32(i/+l)(i/+2) 32(i/-l)(i/-2) l%{v 2 -l) 2 

1 _ q{u+2) q{u-2) q 2 {u+i) , g 2 (t^-4) _ q 2 {p 2 + l)u 

4{u+l) 4(^-1) 32(i/+l)(i/+2) i ~ 32(i/-l)(i/-2) 16(^ 2 — l) 2 

+ 0((7 3 ), (B.7) 



which is called the Mathieu sine function. From this solution we can similarly estimate the 
amplitude of the particle's movement as 

max(|S|) = I (l + ^ + + 0(fl ») , (B.8) 

a = 0,/?= ( ^g +5) , for,<l, 
« = l,/3= 17+ y + 3 2 f + - 3 , forl<^<2, 
a = 1 , p = -i8-8^+2^ + ^ t for v > 2 



For a = 0, it can be expressed as 

8^ 16v^ 



max(|S|) = ^ + - + 0(q 3 ) , (B.9) 



while for a > we obtain 

max(|S|) = -L + r( aq r U + ^2 3/2 + °(9 3 ) . ( B -!0) 
Va 2va(va - 1) 8(a — l) 2 a 6 ' 2 

a = 0,(3= 2(-l)^(a-4)(a-l)+a(a-2) (a+ 5) ? for a < 1 , 

4(^l)lv^](^2-^+2a+a3/2) +a ( 17+ 15^+3a+a3/ 2 ) 

1 _ 2(-l)l^ r l(a-4)(a-l)+a(-18-8 v ^+5a+2a 3 / 2 +a 2 ) 

a = 1 , () = ^ '- , for a < 4 , 

where [\/a\ denotes the integer part of \/a. The above amplitudes tell us how far the 
particle will go, i.e. the absolute value of the largest value that the coordinate \x(t)\ will 
take on, with the initial condition that the particle starts from x(0) = with velocity 
x(0) = 1 at time r = 0. If one wishes to start with the velocity vq, it suffices to multiply 
the solution (B.7) by vq since the equation of motion is linear and hence also the amplitude 
scales accordingly. 
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C A diagrammatic derivation of the effective potential 

Let us consider the classical mechanics problem 

L = ^mX 2 - gX n cos(ut) , (C.l) 

with n > 2. Carrying out a Fourier transform of the field X(t) = j dvX(v)e lvt and 
integrating the Lagrangian over time to obtain the action, we get 

i r . .2 



S = -mjdu v A X(v)X*(v) (C.2) 
- | J dv\dv 2 ■ ■ ■ dv n X(y\)X(y 2 ) ■ ■ ■ X(v n )5(vi + v 2 -\ 1- v n - lj) 

~ 2 J dUldU2 '" dUn X ( u l) X ( u 2) ■ ■ ■ X(v n )5(vi +V 2 ~\ h^+w). 

To obtain the effective action, we will split the field X{v) into two fields according to 
whether \v\ < A <C uj in which case we will call the field still X(v) or whether \v\ > A 
where we call it £(V). This splits the field into the drift part X and the rapidly moving 
part £ and leads to the action 



S=-m[ dv v 2 X(y)X*(y) + -m [ dv v 2 i{v)C{y) (C.3) 

2 J\u\<A 2 J\u\>A 

- ^ / dv\dv 2 ■ ■ ■ dv n X(vi)X(b> 2 ) ■ ■ ■ X{v n )5{ui + v 2 H h v n — uj) 

1 J\vi\<A 

- ^ / dv\dv 2 ■ ■ ■ dv n X(vi)X(b> 2 ) ■ ■ ■ X{v n )5{ui + v 2 H h v n + uj) 

1 J\vi\<A 

- -y \ Ul \ >A dvidv 2 ■■■dv n ^{v{)X{v 2 ) ■ ■ ■ X{u n )5{vi + z/ 2 H \-u n -ui) 

QTi f 

- -y L X | >A dvxdv 2 ■ ■ ■ du n ^(v{)X{v 2 ) ■ ■ • X(z/„)<5(z;i + v 2 H hf„ + w) 

Li|>A dv 1 dv 2 ■■■dv n £{i>i)i{v 2 )X(v?,) ■ ■ ■ X{v n )5{vi + v 2 H h^„-w) 

4 J\vi\>A 
N<A 

j Li\>A dv x dv 2 ■■■dv n £,{vi)£,{v 2 )X{uz) ■ ■ ■ X{y n )8(ui + v 2 -\ h p n + uj) 



N<A 



| / dv x dv 2 ■■■dv n i{yx)i{y 2 ) ■ ■ ■ £{v n )5(vi + u 2 ^ \-v n -ui) 

1 JWi\>A 

| / dv x dv 2 ■■■dv n i{vi)i{h"2) ■ ■ ■ COn)^ + v 2 H \-u n + uj) 

1 J\vA>A 



'H>A 

To get the effective action one solves £ in terms of X and substitutes back into the action. 
Typically the energy v of ^{y) is of order ui so £ is of order 1/uj 2 . The solution of £ in terms 
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of X can then be done iteratively and be presented by diagrams. From this action we can 
read off the propagator as well as the vertices for calculating Feynman diagrams in order 
to construct the effective action in terms of the field X(v). Let us note that the second and 
third lines in the above action tell us that to zeroth order in £ there is no interaction among 
the drift part of the field X as the sum of v\ + u 2 + • • • + v n is also much smaller than u 
unless n is a huge integer. At this point we have established for finite (and reasonable) n, 
the effective action has no potential to order g. Let us now calculate the leading correction 
(at order g 2 ) to the effective potential case by case. 

C.l Quadratic case 

For n = 2 there is a subtlety that there is some sort of mass term. We will not sum up the 
mass insertions but instead consider them as interactions and hence the propagator of £ is 

D = --^, (C.4) 

while the leading vertices are 

© = -igb(i>\ + v 2 + oj) , = -igb{y\ + v 2 - . (C.5) 
To order g 2 there are two tree-level diagrams 

e © © e 

+ 



X(v) X{v) X(v) £(v + uj) X(v) 

which gives the amplitude when taking into account the appropriate symmetry factor 

1 3v 2 



9 2 v 2i \ , 9 2 v2/ \ 9 2 



0J* OJ^ 



X 2 {y). (C.6) 



2m(u! — v) 2 2m(u! + v) 2 m 

To obtain the effective potential we need to make the inverse Fourier transform, such that 
vX becomes X. Finally we get 

VcS = -^X 2 + ^X 2 , (C.7) 
which equals that of eq. (2.11) when identifying g = A 3 . 
C.2 Quartic case 

The n = 4 case is analogous but involves slightly more combinatorics. The propagator of 
£ is still that of eq. (C.4) and the leading vertices are 

© = -H2g5{vi + v 2 + v z + v A + uj) , © = -iUgSfa + v 2 + v z + v A - J) . (C.8) 

To order g 2 there are two tree-level diagrams: 

x{ Vl ) iM) 




£,{V\ + V2 + V 3 - Uj) 
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and the same with © and © interchanged. Taking into account the appropriate symmetry 
factor, the sum of the above diagrams gives the amplitude 



2g 2 



+ 



2g 2 



m(u) — v\ — vi — m(oj + v\ + v 2 + vj) 2 



1 



1 

v 2 



m 



_j_ 3(t/i + v 2 

CO 2 



+ 0(u 



xwxtom^xMmvtixtewv! 



(C.9) 



As before, to obtain the effective potential we need to make the inverse Fourier transform, 
such that vX becomes X. Finally we get 



V, 



off 



raw' 



■X b + 



W8g< 



muj H 



X 4 X' 



(CIO) 



which equals that of (2.14) when identifying g = A 5 . 

The above effective potentials are conservative and independent of time even though 
the original Lagrangian was time-dependent. This effective conservation of energy can be 
understood as follows. By means of our approximation, we have < w. If we look at 
the above tree-level diagram we can see that by starting with a © vertex, the only way 
to obey the constraint that be <C u) is to end the diagram with a © vertex. We can 
construct bigger tree diagrams (to higher order in g and in 1/uj) by inserting for instance a 
vertex with two Xs and two £s (which comes with a g). As long as the number of inserted 
vertices k is small enough we need the same number of © as © vertices 





Hence, to this order the effective action conserves the energy of the drift degrees of freedom 
X and does not have explicit time dependence. This diagram is just an example, many 
other tree-level diagrams of different shape can be constructed. We can now estimate for 
how large k we will be able to end the diagram with © vertex and hence introduce u 
dependence in the effective action. Generically, when the number of vertices is different 
from the number of © vertices, the effective action re-acquires time dependence. If w/A is 
of order of the number of external legs L = 2(3+k) which corresponds to order n; _2 ( fc+1 ) and 
g 2+k , that is to order to/ A — 4 in 1/w, 3 then the effective potential will not be conservative. 
The argument holds for any diagram and not only the presented example. 
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